Can neutralinos in the MSSM and NMSSM scenarios still be light? 
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Since the recent results of direct detection experiments at low mass, many authors have revisited the case 
of light (1 -10) GeV WIMPs. In particular, there have been a few attempts to explain the results from the 
DAMA/LIBRA, CDMS and/or CoGeNT experiments by invoking neutralinos lighter than 15 GeV. Here we 
show that in the MSSM, such light particles are completely ruled out by the TEVATRON limits on the mass 
of the pseudoscalar Higgs. On the contrary, in the NMSSM, we find that light neutralinos could still be viable 
candidates. In fact, in some cases, they may even have an elastic scattering cross section on nucleons in the 
range that is needed to explain either the DAMA/LIBRA, CoGeNT or CDMS recent results. Finally, we revisit 
the lowest limit on the neutralino mass in the MSSM and find that neutralinos should be heavier than ^ 28 GeV 
to evade present experimental bounds. 



I. INTRODUCTION 

Recently, both the CDMS and CoGeNT experiments have an- 
nounced results which have been regarded as possible hints of 
the existence of dark matter particles. In particular, the CDMS 
experiment has claimed the detection of two "anomalous" 
events in a blind analysis 1 1 1 while the CoGeNT experiment 
has observed an unexplained rise in a p-type point contact 
(PPC) spectrum devoid of most surface events O. Together, 
these results might confirm the long-standing claims of detec- 
tion of dark matter annual modulation by the DAMA/LIBRA 
experiment. Therefore, these findings deserve some attention. 

The most puzzling aspect of these recent claims is that they 
seem compatible with the existence of rather light dark matter 
particles. For example, the assumption of dark matter par- 
ticles with a mass between 7 and 1 1 GeV seems to provide a 
good fit to CoGeNT data (even though the null hypothesis also 
provides a similar fit) while particles from 7 GeV to 40 GeV 
could explain the CDMS events. This mass range is also com- 
patible with the DAMA/LIBRA O Hi findings which favour 
particles in the 10-15 GeV range assuming spin-independent 
(SI) interactions and no channeling. 

Although these claims sound quite exciting, they raise several 
issues. First of all, CDMS claimed detection of two events 
only. This is not statistically significant to enable any firm 
conclusion related to dark matter. However, detection of more 
than five events would have been in conflict with the nega- 
tive results of XENON 10 [5J. Secondly, the possible inter- 
pretation of a ~ 7 — 1 1 GeV dark matter candidate with a SI 
cross section around 10""^^ — 10~'^^cm^(as proposed by the 
CoGeNT collaboration) is severely challenged by the negative 
results of the XENON 100 experiment |6|. Nevertheless, it 
was also demonstrated that the uncertainties on the dark mat- 
ter escape velocity Vesc and scintillation efficiency Leff may 
offer a way to reconcile both XENON 100 and CoGeNT re- 
sults HHHl. Finally, the CDMS and CoGeNT events seem to 
exclude the high value of the spin-(in)dependent cross section 
that is favoured by the DAMA/LIBRA experiment. 

Nevertheless, it is worth mentioning that the mass range that 



is under consideration is not excluded neither by the ISR con- 
straints at LEP ifTSl nor by the potential Sunyaev-Zel'dovich 
effect that might be generated by relatively light dark mat- 
ter particles in clusters of galaxies. (161 [TTl. Although it 
is premature to draw any conclusion about the nature of the 
DAMA/LIBRA, CDMS and CoGeNT results, it is certainly 
worth reinvestigating the possible dark matter candidates that 
are expected in the low mass region. 

In this paper, we shall only consider supersymmetric candi- 
dates and in fact focus more specifically on the neutralino 
(Xo)- We are seeking a light Xo with a relatively high elas- 
tic scattering cross-section with matter. In order to retain the 
connection with the "dark matter" interpretation, one should 
also require that this candidate has a relic density which rep- 
resents a subsequent fraction (if not all) of the dark matter 
abundance inferred from the combination of CMB data with 
other cosmological (e.g. SN, BAO) observations 1 18|. 
Recently, the authors of Ref. 1 19 1 claimed that light neutrali- 
nos in MSSM-EWSB scenarios with non-universal gaugino 
masses could actually do the job. This work was followed 
by other claims (using supersymmetric extensions of the 
Standard Model ll2Q]425l ) to explain CoGeNT and/or CDMS 
events. However, light MSSM neutralinos were in fact inves- 
tigated before these experimental claims |T5lfT9, 26 - 28] and 
just after the study of the astrophysical signatures that are to 
be expected from light WIMPs |29|. Neutralino masses down 
to ~ 5 GeV were found. However in Ref. (301, it was pointed 
out that the improved measurement of Bs jn/n excludes such 
low neutralino masses. Very light neutralinos that would con- 
stitute hot dark matter were however found to be consistent 
with all experimental constraints |31 1. 

Here, we reinvestigate the MSSM in light of the latest TEVA- 
TRON results to determine whether light neutralinos can ex- 
plain the direct detection signals or not. We find that, when 
taking into account the Tevatron results on the supersymmet- 
ric Higgs at large values of tan p, the lower limit on the mass 
of the neutralino increases to 28 GeV. We then consider an ex- 
tension of the MSSM with an extra singlet, the NMSSM. Neu- 
tralino dark matter in singlet extensions of the MSSM were 
first investigated in Ref. |32|. When the singlet Higgses are 
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light we find scenarios with neutralinos of a few GeV that 
evade all constraints and yet predict a direct detection signal 
in the region preferred by recent experiments. 
In section II, we describe the method used for exploring the 
parameter space. We analyse scenarios with light neutralinos 
in the MSSM assuming different priors in section III while the 
results for the NMSSM are presented in section IV. 

II. METHOD 

To efficiently explore the multi-dimensional parameter space, 
we have performed a Markov Chain Monte Carlo analysis 
(based on the Metropolis-Hastings algorithm). We have used 
micrOMEGAs2.4 |[33l |34| to compute all observables. This 
code relies in turn on SuSpect |35 1 for calculating the par- 
ticle spectrum in the MSSM and on NMSSMTools |[36]| for 
calculating the particle spectrum and the various collider and 
B-physics constraints in the NMSSM. 

We use the method of burn-in chains, i.e. we first explore the 
parameter space till we find a point with a non- vanishing like- 
lihood. When such a point is found, we continue the chains, 
keeping all the points that are retained by the MCMC. How- 
ever, since it is difficult and time-consuming to find a good 
starting point, we require to speed up the process that the like- 
lihood times the prior (hereafter referred to as Q) associated 
with the starting point exceeds the value Q> 10"^^, and use 
an exponential prior on m^^ to make sure that the starting point 
is within close proximity of the low neutralino mass region. 
However, when this point is found, we replace the exponen- 
tial prior on m^^ by a flat prior. Since low mass neutralinos 
are quite unlikely with respect to heavier ones (and since find- 
ing them also requires a certain amount of fine-tuning), we 
have decided to perform two independent scans. One aims 
at exploring the mass region ranging from to 15 GeV (this 
range of mass is particularly relevant for the CoGeNT and 
DAM A/LIBRA interpretation) and the second one aims at 
exploring the range from to 50 GeV (to include the pre- 
ferred region of the two WIMP recoil-like events reported by 
CDMS-II). A proper exploration of the parameter space is ob- 
tained after generating approximately 50 chains of 10^ points 
each. 

The total likelihood function for each point is the product of 
the likelihood functions evaluating the goodness-of-fit to all 
the data set that are displayed in Table. |l] These include B 
physics observables, the anomalous magnetic moment of the 
muon, {g — 2)^, the Higgs and sparticles masses obtained from 
LEP and the corrections to the p parameter. For the MSSM 
case, only LEP mass limits on new particles were taken as a 
sharp discriminating criterion with X = or 1 . Other crite- 
ria had some tolerance. For the NMSSM, limits on the Higgs 
sector, on the Z partial width and on neutralino production as 
computed by NMSSMTools were also taken as a sharp dis- 
criminating criterion. 

We use a Gaussian distribution for all observables with a pre- 
ferred value ju±a, 

F2{x,lA,G)=e 2a2 (1) 



and 

^3(x,A/,a) = ^z^. (2) 

for observables which only have lower or upper bounds. The 
tolerance, a, is negative (positive) when one deals with an 
upper (lower) bound. 



TABLE I: List of constraints, from Ref. |[37l unless noted otherwise 



constraint 


value/range 


tolerance 


applied 


Smasses 




none 


both 


^WMAph^ 


0.01131-0.1131 


0.0034 


both 




25.5 10-10 


Stat: 6.3 lO-^^ 
sys: 4.9 10-^^ 


both 


Ap 


< 0.002 


0.0001 


MSSM 




3.52 10-4 [g8l[39l 


th: 0.24 10-4 
exp: 0.23 10-^ 


both 




<4.7 10-^ 


4.7 10-10 


both 


R{B TV) 


1.28 |38 | 


0.38 


both 


rriH 


> 114.4 


1% 


MSSM 




< 1.7 MeV 


0.3 MeV 
none 


MSSM 
NMSSM 


e^e- XiX2,3 


<0.1pb |40| 


0.001 pb 
none 


MSSM 
NMSSM 


AM, 


117.0 10-1^ GeV 


th: 21.1 10-1^ GeV 
exp: 0.8 10-^^ GeV 


NMSSM 


AM^ 


3.337 10-1^ GeV 


th: 1.251 10-1^ GeV 
exp: 0.033 10- GeV 


NMSSM 



Finally we also require that the neutralino relic density satis- 
fies 

100% Q^wMAph^ > ^jh^ > 10% awMAph^^ (3) 

with Q^wMAph^ = 0.1131 ±0.0034 (HI. The cases where 
Q,^ < Q^wMAP should correspond to scenarios in which there 
is either another (if not several) type of dark matter particles 
in the galactic halo B2ll or a modification of gravity (cf e.g. 
1 43 1). In case of a multi component dark matter scenario, there 
could be either very light e.g. 1 44- 47 1 or very heavy particles 
(including very heavy neutralinos), depending on the findings 
of direct detection experiments. 



III. MSSM SCENARIOS 

In what follows, we consider the MSSM with input parame- 
ters defined at the weak scale. We assume minimal flavour 
violation and equality of the soft masses between sfermion 
generations. We further assume a common mass mj for all 
sleptons, and for all squarks niq (but we have checked that we 
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found consistent results by relaxing this universality assump- 
tion). We allow for only one non-zero trilinear coupling, At. 
The gaugino masses Mi and M2 are free parameters which, 
in particular, allows to have Mi <C M2 implying a light neu- 
tralino much below the EW scale. The parameter M3 satisfies 
the usual universality condition in GUT scale model, that is 
M3 = 3M2. The Higgs bilinear term, ju, the ratio of Higgs 
vev's, tanP and the pseudoscalar mass Ma are also free pa- 
rameters. 

This MSSM-EWSB model with only eight parameters can re- 
produce the salient features of neutralino dark matter. Indeed, 
apart from the mass of the LSP, the most important parame- 
ters are the gaugino/higgsino content of the LSP, determined 
by jn and Mi , M2, tan p, as well as the mass of the pseudoscalar 
which can enhance significantly neutralino annihilations into 
fermion pairs. 



A. Neutralino masses smaller than 15 GeV 

To sample the low neutralino mass range, we take our priors 
in the range 

Ml e [1, 100]GeV M2 e [100,2000]GeV 

jn e [0.5, 1000]GeV tanP e [1,75] 

mj e [100, 2000]GeV G [300, 2000]GeV 

At e [-3000, 3000] GeV ma e [100, 1000]GeV 

We consider separately the cases ^ > and ju <0. The results 
of our MCMC simulations for ^ > are displayed in Fig. [T] 
and [21 

Fig. [T] represents Q/Qmax^ the weigth normalized to the best 
weight, with respect to the free parameters that we have con- 
sidered. The first plot shows that the bino mass is peaked 
around Mi G [15,19] GeV while the second plot shows that ju 
is below 150 GeV. That is, it is near the lower bound that satis- 
fies the LEP limits on charginos. Thus, the LSP is dominantly 
bino-like with a small Higgsino component. 
The third and fourth plots in Fig. [T] show that tanP is very 
large (tanP G [40,60]) and ma is relatively small (rriA G 
[120, 170] GeV). This basically indicates that the main neu- 
tralino pair annihilation proceeds through the s -channel ex- 
change of a light pseudo-scalar Higgs boson. The results 
also show that the sleptons and squarks are preferably heavy 
(mf G [500, 1200] GeV and G [0.8,2] TeV). 
In Fig. [2] we display the same quantity but with respect to the 
neutralino mass. As one can see, the preferred value for the 
neutralino mass lies in between 13 and 15 GeV. We found 
no scenario where the neutralino mass would be smaller than 
10 GeV. 

In Fig. [3] we display the (spin-independent) elastic scattering 
cross-section a^^ times the fraction of neutralino in the halo 
^ versus the neutralino mass and the limits from CDMS and 
XENON 100. Here (and in the following), we have assumed 
values for the quarks coefficients in the nucleon (defined by 
setting CnN = 45MeV, Cq = 40MeV in micrOMEGAs |[48l ) 
that lead to rather low cross-sections in order to be conserva- 



tive in our predictions . Since there are uncertainties on the 
escape velocity and scintillation function of XENON 100, we 
also performed a rescaling of L^/f with the energy (see B9ll ) 
and kept a conservative energy-dependent value for Le/f. In 
principle, this should enable us to derive the most conservative 
limit as possible. 

Since all the light candidates that we have found lie within 
the excluded region (which is defined by both CDMS and 
XENON 100 exclusion curves), our results strongly suggest 
light MSSM neutralinos cannot explain neither the CoGeNT 
nor the DAMA/LIBRA data. Hence, if dark matter is indeed 
as light as 10 GeV, MSSM neutralinos are extremely unlikely 
to be a viable explanation. 

We reach the same conclusion by using collider constraints. 
Indeed, since the points that we have found are associated with 
a light pseudo-scalar Higgs, one can use the 95% exclusion 
limits from TEVATRON in the MSSM plane of tanP versus 
150] to constrain the scenarios selected by our MCMC. 
These collider limits were obtained by studying the MSSM 
Higgs boson production in association with b quarks in the 
x+x~ final states with up to 2.2 fb~^ of data. Two benchmark 
scenarios were in fact considered to describe the situations of 
maximum and no mixing in the stop sector 15TI for both /n>0 
and /i < 0. 

In Fig. [4] we thus plot the results of our MCMC analysis and 
superimpose the TEVATRON constraints. We find that, in- 
dependently of the direct detection constraints, none of the 
scenarios corresponding to light MSSM neutralinos and ^ > 
can survive the TEVATRON constraints on the mass of the 
pseudo-scalar Higgs. The fact that neutralinos lighter than 15 
GeV (/i > 0) are excluded by both direct detection and collider 
constraints enable us to conclude that MSSM-EWSB neutrali- 
nos (in scenarios with jn > 0) must be heavier than 15 GeV. 
This argument is actually valid whatever the precise value of 
the relic density. This is therefore a very strong conclusion. 
Note that our results are somewhat discrepant with |[T9ll . This 
may originate from the fact that Ref. |19| did not take into 
account the latest constraint on Bs /njn. Indeed, we do obtain 
many more points at low mass if we disregard this constraint, 
as was also pointed out in Ref. |30|. However, all these "new" 
points are also excluded by the TEVATRON limits. 
We have performed a similar analysis for the case ^ < 0. How- 
ever, we found that it is extremely difficult to obtain a good 
starting point. Besides the total Likelihood of the points re- 
tained by the MCMC is much smaller than that for /i > 0, 
this is mainly due to the {g — 2)^ constraint and to a lesser ex- 
tend to the constraints from B -physics. Hence, MSSM-EWSB 
neutralinos lighter than 15 GeV are also ruled out in the case 
A/<0. 2 

Of course our lower bound on Q^^^^ (see Eq 31) is arbitrary. 

ofligl] 



One could consider an even smaller fraction < 



ight neutrali- 



^ The elastic scattering cross-section can be up to one order of magnitude 

larger for other choices of the quark coefficients. 
^ In Ref. 1 52 1 it was shown that the {g — 2)^ constraint could be avoided if 

one takes opposite signs for gaugino masses with both jj <0 and M2 < 0. 

We have not considered this class of scenarios. 
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FIG. 1: MSSM-EWSB scenario for > and < 15 GeV. These plots represent the rescaled weight Q/Qmax of the points selected by the 
MCMC versus the free parameters that we have considered. Curves in dark blue correspond to points with a likelihood greater than 99.4 %; 
Curves in blue, correspond to points with likelihood greater than 95.4 % and smaller than 99.4 % of the maximum Likelihood and points in 
pale blue are all the remaining points having a likelihood greater than 68 %. 
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FIG. 2: Q/Qmax with respect to the neutralino mass in MSSM- 
EWSB scenario for > and m^^ < 15 GeV. We use the same color 
code as in Fig.[T] 



nos in the halo. In this case, not only would they be a minor 
source of dark matter but also their energy density distribu- 
tion in the halo may be affected by the nature of the main 
candidate; this is hard to estimate. We have checked never- 
theless that taking Q.^ <C 10%Q.wmap does not change our re- 
sults. This can be understood easily. Taking Q,^ <C Q^wmap 
requires a large pair annihilation cross-section (cJv^^^^). Yet, 
GVj^^ is proportional to the neutralino mass squared. Hence, 
the lightest the neutralino, the smallest the pair annihilation 
cross-section, that is the largest value of Q.'^h^. Thus, one can- 
not take the neutralino mass arbitrarily small. To enhance the 
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m^o [GeV] 

FIG. 3: MSSM-EWSB scenario with > and < 15 GeV. Spin- 
independent cross section on proton times the fraction of neutrali- 
nos in the Milky Way dark halo (^) versus the neutralino mass m^. 
The dark red (light pink) points have a likelihood greater than 99.4% 
(68%). 



cross-section, the sole viable option is to invoke neutralino 
pair annihilations through the exchange of a relatively light 
pseudo-scalar Higgs (which is in agreement with our findings) 
since neutralino co-annihilations with the Next-to-LSP is im- 
possible owing to the smallness of the neutralino mass. 
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FIG. 4: Distribution of the points selected by our MCMC analysis in 
the tan P — plane in the MSSM-EWSB scenario with ^ > and 
rur^ < \5 GeV. The TEVATRON limits are displayed for the case of 
no-mixing (dash) or maximum mixing (full) in the stop sector, same 
color code as in Fig.|3] 



B. Neutralino masses less than 50 GeV 



towards higher values. Furthermore tanP varies in a wider 
range since the 5-physics constraints are relaxed at large val- 
ues of niA. 

In Fig. [6j we display the points selected by the MCMC in 
the plane (mA, tan P) where we have superimposed the TEVA- 
TRON constraints. Interestingly enough, at low values of the 
mass of the pseudo-scalar Higgs rriA < 300 GeV, there are two 
separate regions of tan p. One is peaked around 10-20 while 
the second lies between 50 — 70. Typically when the pseu- 
doscalar is light, constraints on 5-physics decrease the value 
of the likelihood especially when tanP is large. However as 
we have seen previously the new channel for neutralino an- 
nihilation through a pseudoscalar exchange leads to an ac- 
ceptable relic density and to a good global likelihood when 
tanP > 50. Note that, even though very large values of tanP 
do not appear plausible, they do indicate the type of regions 
that lead to a neutralino mass in the [1 , 50] GeV range. 
The spin-independent cross section versus the neutralino mass 
is displayed in Fig. [7] The scenarios where the neutralino 
mass is greater than 28 GeV survive both the TEVATRON 
and Direct Detection limits. Although such a value is likely to 
be irrelevant to explain CoGeNT data, it might be important 
in light of the two CDMS "events". 



Given our conclusion concerning light neutralinos, it is inter- 
esting to derive the lower limit for the neutralino mass in the 
MSSM. For this, we focus on scenarios where the neutralino 
mass ranges from 1 to 50 GeV. 




FIG. 5: Same as in Fig.[T]for the MSSM-EWSB scenario with in>0 
and < 50 GeV 



The distribution of Q/Qmax for the free parameters of the 
model are summarized in Fig|5] As one can see, the preferred 
value of the neutralino mass lies in the [28,50] GeV range 
and is similar to the value of M\ which lies in [40, 50] GeV. 
Since M\ < ju, such neutralinos are mostly bino-like. Never- 
theless some Higgsino component is necessary for efficient 
neutralino annihilation through Z or light Higgs exchange, 
hence the preference for small values of jn. As compared to 
the previous case, the pseudoscalar Higgs is no longer neces- 
sary to provide efficient annihilation, and therefore ma moves 
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FIG. 6: Distribution of the points selected by our MCMC analysis 
in the tanP — ma plane in the MSSM-EWSB scenario with ju > 
and < 50 GeV. In red, we display the points which are excluded 
by both TEVATRON, XENON 100 and CDMS. In yellow, we show 
the points which satisfy TEVATRON and which are excluded by 
XENON 100 and CDMS and in green, all the points which survive 
both constraints. 

For /i < 0, we again found that the Likelihood of the points 
is much smaller than for /i > 0, owing mainly to the — 2)^ 
constraint which disfavours large values of tan |3 as well as low 
values of M2 and mj 1 15 1. 



IV. NMSSM SCENARIOS 

The Next-to-Minimal Supersymmetric Standard Model 
(NMSSM) is a simple extension of the MSSM that provides 
a solution to the naturalness problem. This is achieved by 
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FIG. 7: Spin independent cross section vs the neutralino mass in the 
MSSM-EWSB scenario with ^ > and m^^ < 50 GeV, same color 
code as in Fig. [6] 



the introduction of a gauge singlet superfield, denoted by S. 
The VEV of this singlet determines the effective parameter 
jLi = ^{S) which is then naturally of the EW scale |53|. The 
part of the superpotential involving Higgs fiels reads 



W 



(4) 



and the soft Lagrangian 



(5) 



The NMSSM contains three neutral scalar fields, h\ , /i2, ^3 and 
two pseudoscalar neutral fields, a\^a2 as well as a charged 
Higgs, H^. The model also contains five neutralinos, the new 
field is the singlino, S. For a pure state the singlino mass is 
simply 1.5 3 J 



(6) 



After using the minimization conditions of the Higgs poten- 
tial, the Higgs sector is described by six free parameters, 
^,tanP as well as X^k^A^^A^. Other free parameters of the 
model are, as in the MSSM, the soft masses for sfermions, 
trilinear couplings and gaugino masses. 

An important feature of the model is that both the singlino 
and the singlet fields can be very light and yet escape the LEP 
bounds. This is because these fields mostly decouple from 
the SM fields 1531 . This opens up the possibility for new an- 
nihilation mechanisms for light neutralinos in particular if the 
LSP possesses an important singlino component. The singlino 
can annihilate efficiently through the exchange of light singlet 
Higgses as well as into light Higgs singlets L54J . 



A. Neutralino masses smaller than 15 GeV 

To explore the parameter space of the NMSSM model that 
allow for light neutralinos we follow the same procedure as 
for the MSSM. Our priors lie in the range: 



Ml G [1,200] GeV 
A/G [0.,1000]GeV 
Xe [0,0.75] 



M2 e [100,2000]GeV 
tanPG [0.1,65] 
KG [0.,0.65] 



Ax G [-2000, 5000] GeV A^ G [-5000, 2000] GeV 

mj G [100, 2000]GeV G [300, 2000]GeV 

At G [-3000, 3000] GeV 

As before, we assume common soft masses for squarks and 
sleptons and we keep the gaugino masses Mi and M2 uncor- 
rected while M3 = 3M2. 

The distribution of Q/Qmax for the free parameters are dis- 
played in Fig. [8] As in the MSSM the bino mass is peaked 
below 20GeV although a long tail extends to 200GeV. In this 
tail the LSP is mostly a singlino which mass is determined 
from Eq. [6] The parameter X that determines the mixing of 
the singlino to other neutralinos is never very small, so that the 
singlino does not decouple completely. The preferred values 
for /J ^ 150- 250 GeV are higher than in the MSSM. On the 
one hand LEP2 limits on e^e~ x^x^ or on the light Higgs 
constrain low values of ju while a light singlino LSP prefers 
low values for jn, Eq.|6] The parameter K << 1 also favours a 
light singlino. Intermediate values of tan P are preferred. The 
parameter Ak that controls the mass of the singlet Higgses is 
always small to ensure a light scalar/pseudoscalar as required 
for LSP annihilation while A^ is usually well above 1 TeV. 
Sleptons are preferably light while squarks are above ITeV. 
The LSP mass ranges from l-15GeV with a distribution 
peaked towards higher masses. This LSP is either mostly bino 
or mostly singlino with in any case some higgsino component. 
The most important feature of this scenario is the fact that 
the Higgs spectrum is constrained: one always predict a light 
scalar, dominantly singlet, with a mass below 120 GeV (gen- 
erally below 30GeV) as well as a pseudoscalar singlet with 
a mass preferably below 30GeV, see fig. fTOl Note that the 



value of 30GeV for the mass corresponds to twice the neu- 
tralino mass and is thus just a consequence of the prior on 
the neutralino mass. Furthermore we find generally that either 



,/2<l. 



1, cf Fig. 11 ) or that m^o > ruh^ 



4 GeV (with a similar mass splitting with 
This is because the anni- 



hilation of lEe light LSP relies either on pseudoscalar/scalar 
exchange or on the new light scalar pairs final states. The rest 
of the Higgs sector consists of MSSM-like doublets with pre- 
ferred values for the heavy neutral and charged scalars above 
2TeV. Note that we have checked that the recent re-analysis of 
LEP2 limits on a Higgs decaying into two light pseudo- scalars 
did not put further constraints on our model parameters (551 . 
The light LSP scenarios can be classified in three broad 
classes: 1) a (pure or mixed) singlino LSP annihilating via 
pseudoscalar/scalar singlet Higgses into fermion final states, 
for this only a small singlino component of the LSP is neces- 
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FIG. 8: Q/Qmax for the points selected by the MCMC versus the value of the free parameters that we have considered in the NMSSM for > 
and <\5 GeV, same color code as Fig.l 



sary. 2) a bino LSP with small higgsino/singlino components 
annihilating into a pair of light scalar Higgses or 3) as in the 
MSSM a bino LSP with some Higgsino component annihilat- 
ing via Higgs doublets. This channel is more efficient at large 
values of tanP although the B-physics constraints severely re- 
strict the very large values of tanp. 

The predictions for the elastic scattering cross section span 



several orders of magnitude, from 10 to 10 



see 



Fig. [12] The largest cross sections are found in scenarios with 
a light hi, for example for a^^ > 10~^^(10~^^)cm^ requires 
< 20(8) GeV. ^ At first sight this can be a bit surprising 



m 

since such a light Higgs is dominantly singlet and thus couples 
very weakly to quarks in the nucleon - recall that the hiqq cou- 
pling is only possible through the doublet component- nev- 
ertheless this suppressed coupling is compensated by an en- 
hancement factor due to the small hi mass in the propaga- 
tor, oc l/ml^. In scenarios where the elastic scattering cross- 
section is large, the LSP is generally dominantly bino (or, in 
a few cases, a singlino) with a non-negligible higgsino frac- 
tion. This means that the doublet /z2 also contributes to the 
spin independent cross section since the LSP coupling to the 



doublet depends on the Higgsino component of the LSP. ^ The 
very low cross-sections are found in scenarios where the LSP 
pair-annihilation benefits from the enhancement of the pseu- 
doscalar exchange in the s-channel near the resonance while 
the elastic scattering cross-section, which proceeds through 
scalar exchange in t-channel, does not benefit from a similar 
enhancement. 



^ Note that in Ref (56| an upper bound on the SI cross-section was obtained 
by optimizing the contribution of the doublet exchange, we found larger 
cross sections because we explored regions where the additional contribu- 
tion of the light singlet was important. 



^ In Ref. 1 22 1, light CP-even Higgs exchange was proposed to explain the 
CoGeNT result with light neutralinos within the BMSSM model. 
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FIG. 9: Q/Qmax with respect to the neutraUno mass in the NMSSM 
for ^ > and m^^ < 15 GeV, same color code as in Fig.fl] 




[GeV] U^^ [GeV] 

FIG. 10: Q/Qmax with respect to the mass of the scalar (left) and 
pseudoscalar Higgs (right) in the NMSSM for ^ > and m^^ < 15 
GeV, same color code as in Fig.[T] 




m^o [GeV] 

FIG. 12: NMSSM scenario with > and m^^ < 15 GeV. Spin- 
independent neutralino-proton cross section times the fraction of 
neutralinos in the Millcy Way darlc halo (^) versus the neutrahno mass 
rrij^. We use the same color code as in Fig.fl] The CDMS limits cor- 
respond to the plain curve while the XENON limits correspond to 
the dash curve. 



low the present limits. Note that we also found that neutrali- 
nos lighter than 15 GeV were accompanied by a singlet pseu- 
doscalar having a mass ^ 2m^. However, these scenarios did 
not have a significant singlino component. In addition, they 
had a lower likelihood than in the case of heavier neutralinos. 



6000 



4000 



2000 




-20 
Mn-Mn 12 [GeV] 

1 



FIG. 1 1 : NMSSM scenario for ;t/ > and < 15 GeV. Mass differ- 
ence rrij^ — Mai We use the same color code as in Fig.[T] 

We have also explored the region of parameter space which 
gives neutralino masses up to 50GeV. As for the MSSM we 
found many scenarios that satisfies all the constraints. These 
have a neutralino mass near Mz / 2 and a SI cross section be- 



B. LHC signatures for light neutralino scenarios 



Light NMSSM neutralino scenarios (m^ < 15 GeV) can lead 
to distinctive signatures at colliders both in the Higgs and the 
neutralino sector. First, the usual dominant decay mode of the 
light scalar doublet Higgs can be greatly suppressed because 
of new decay modes into neutralinos h2 X^X^ as well as into 
light Higgses h2 h\h\^a\a\. The light scalar/pseudoscalar 
are expected to decay predominantly into bb^ x+x~ or invisibly 
into X^X^. Therefore both the search for the invisible Higgs in 
the W fusion channel with a signature in two tagged jets and 
missing Et and the search for Higgs via Wh2 production with 
h2^h\h\ (571 [58l are important channels at the LHC. 

Distinctive signatures are also expected in the neutralino and 
chargino sectors. The NLSP, X2^ can contain a small compo- 
nent of singlino and therefore be much lighter than Xi- The 
favourite discovery channel, X2 ~^ however much sup- 
pressed since the decay modes can involve light Higgs states, 
X2 ~^ leading to X2 ~^ Xi^^ completely 

invisible decay X2 ~^ XiXiXi- Furthermore the single lepton 
signature of the chargino will be much suppressed since the 
decay W^xH is generally accessible because of the large 
mass splitting with the LSP. 
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V. CONCLUSION 

We have investigated whether Hght neutralinos in the MSSM 
and the NMSSM could still be viable candidates and fur- 
thermore lead to large values of the elastic scattering cross- 
section as required to explain recent CDMS, CoGeNT, 
DAMA/LIBRA data. 

In the MSSM, we basically exclude neutralinos lighter than 
15 GeV. In order to satisfy the WMAP constraint, these parti- 
cles must exchange a light pseudoscalar Higgs at large values 
of tan p. However, recent TEVATRON results on supersym- 
metric Higgs searches rule out this possibility. In addition, 
we found (after imposing all collider and direct detection con- 
straints) that the lower limit on the neutralino mass is about 
28 GeV in the MSSM {jj > 0). 

In the NMSSM, it is easier to find light neutralinos that sat- 
isfy all the constraints. This is generally achieved when 
the LSP contains a singlino component. One salient feature 
is that these neutralinos are accompanied by a light pseu- 
doscalar and/or scalar singlet. While the elastic scattering 
cross-sections are generally much below the reach of cur- 
rent detectors in most scenarios, we found that some points 
can satisfy all the constraints and yet predict elastic scatter- 
ing cross- sections in the "CoGeNT" region. This is achieved 
through the presence of a O(GeV) scalar Higgs. However we 
should stress that the light neutralinos typically do not provide 
as good a fit to the data as the ones around 50 GeV. 



NMSSM scenarios with light neutralinos could be further 
probed with direct detection experiments of increased sensi- 
tivity. They also lead to distinctive signatures in Higgs and 
SUSY particles searches at colliders that could be studied. 
Finally, indirect signatures of the light neutralino (including 
synchrotron emission, see Ref. (591) will be investigated in a 
future work. 

Note added: After submitting this paper, we became aware of 
another preprint 1 60 1 which also finds that light neutralinos in 
the NMSSM can have large elastic scattering cross sections 
through the exchange of a light scalar Higgs. 
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